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SUMMARY 

Based on experimental observations, a mechanism is presented to explain 
the existence of the different morphologies of electrodeposited zinc in alka- 
line solution. The high current density dendrites appear to be due to more 
rapid growth on the non-basal crystallographic planes than on the basal 
plane. The low current density moss apparently results from dissolution from 
the non-basal planes at low cathodic voltages. 

Electrochemical models were sought which would produce such a phenomenon. 
The fundamental plating mechanism alone can account only for different rates 
on different planes, not for zinc dissolution from a plane in the cathodic 
region. Fourteen models were explored; two models were found to be in accord 
with the proposed mechanism. One involves rapid disproportionation of the 
zinc +1 species on the non-basal planes. The other involves a redox reaction 
(corrosion) between the zinc-zincate and hydrogen-water systems. 

INTRODUCTION 

An important aspect of the alkaline zinc electrode is the morphology of 
the zinc which is electrodeposited during charging (refs. 1 and 2). At low 
current density and high zincate concentration, the morphology of the zinc is 
termed mossy. This zinc is in the form of fine whiskers reported to be of 
micrometer diameter (ref. 3). At higher current densities and low zincate 
concentrations, boulders (nodules) of zinc form (ref. 4). These are macro- 
crystalline. At extremely low average current densities, a smooth plate of 
zinc is attainable (ref. 5) and at extremely high current densities, dendritic 
zinc is formed. The morphology is also affected by the presence of impurities 
(ref. 6). 

Because the boulder and dendritic morphologies lead to shorting, the mossy 
form is generally preferred for battery applications (ref. 7). It would be 
very useful to preserve the mossy form at higher current densities in order to 
increase the rate of charging without favoring shorting. A knowledge of the 
mechanism of the mossy morphology would be highly desireable. Generally, 
the morphology is examined from the viewpoint of why dendrites grow (refs. 5 
and 8). 

The purpose of this report is to present and discuss a possible mechanism 
for the formation of the various zinc morphologies and propose several models 
that can give rise to this mechanism. Emphasis is placed on an explanation of 
the mossy deposit. Many of the concepts presented here should be applicable 
to the plating of other metals. 

According to our proposed mechanism, the mossy zinc morphology is due to 
the dissolution of zinc from the non-basal crystallographic planes at low 
cathodic voltages. Two models were found to be in accord with this mech- 
anism. One involves rapid disproportionation of the zinc +1 species on the 






non-basal planes. The other involves a redox reaction between the zinc- 
zincate and hydrogen-water systems. 

OBSERVATIONS LEADING TO THE MECHANISM 

Figure 1 is a Scanning Electron Microscope (SEM) view of dendritic zinc 
deposited at high voltage and high current density. Note that the dendrites 
are composed of shingled hexagonal plates. It seems clear that the faces of 
these plates are the basal planes of the hexagonal zinc crystal lattice. 

The edges of these plates are non-basal planes. One concludes that in the 
dendritic deposition region, zinc plates much more rapidly on the non-basal 
iTOO, lOTO, and OlTo planes than on the basal faces for thin plates to form. 

In this report, reference to non-basal planes and N planes indicate the iTOO, 
lOTO, and OlTO planes; reference to the B plane refers to the basal plane. 

Figure 2(a) is an SEM view of mossy zinc deposited at low voltage and low 
current density. Unlike the dominating two dimensional growth of dendrites, 
the moss grows predominantly in one dimension. The lack of facets or other 
crystal-like appearance suggests that the zinc may not be crystalline. We 
have, however, obtained X-ray and electron diffraction patterns that show the 
moss strands to be the hexagonal zinc crystal structure. Figures 2(b) and (c) 
are respectively a Transmission Electron Microscope (TEM) view of a mossy zinc 
strand and an electron diffraction pattern taken on that strand. The pattern 
in figure 2(c) proves that the strand shown in figure 2(b) is single crystal 
zinc and that the cross section of the strand coincides with the basal plane. 
This orientation indicates that for a mossy deposit, zinc plates much more 
rapidly on the basal plane than on the non-basal sides of the strands. There- 
fore, it appears at first that the relative rates of plating on the basal and 
non-basal planes determine which morphology (dendrites or moss) will form. 

However, two observations concerning the moss seem somewhat contra- 
dictory. They are: (1) the strands of moss are single crystals and (2) the 

SEM pictures never display any typical crystal growth facets. A way to re- 
solve this apparent contradiction is to assume that zinc is, in fact, dis- 
solving off the non-basal planes during the deposition of the mossy zinc: 
dissolving crystals tend to lose their faceting and become more rounded. 

Figure 3 is an SEM picture of zinc crystal growth at conditions inter- 
mediate between the high voltage (non-basal dominating) dendrites, and the low 
voltage (basal dominating) moss. These "boulders" display typical hexagonal 
crystal shape with comparable plating rates on basal and non-basal planes. 

To explain these observations, our proposed mechanism is the following. 

At low cathodic overvoltages, zinc plates on the basal plane but dissolves off 
the non-basal plane, giving rise to the whiskers called moss. At intermediate 
voltages, zinc begins to plate on both sets of planes. The rates on both sets 
of planes are comparable so that boulders result. At high voltages, the rate 
of deposit on the non-basal planes exceed that of the basal plane. The den- 
drites result. 


THEORETICAL MODELS 

The next step is to examine models to find one or more which can give rise 
to the proposed mechanism. The models consist of the electrochemical and 
chemical reactions which are (or may be) involved as well as the applicable 
Butler-Volmer equations; in general two sets are required, one for the basal 
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and one for the non-basal planes. The chief test for a model is its ability 
to predict dissolution of zinc from the non-basal plane at low cathodic 
voltages. 

Fourteen models (A to N) are examined. All of these assume steady state 
conditions, no concentration gradients between the bulk and plating surfaces, 
and a constant ratio between the B and N plane areas. Model A consists of the 
use of only the fundamental equations involved with the deposition of zinc and 
diffusion of species between the N and B planes. Each of the other models 
presented here use the equations operative in model A plus one additional pro- 
posed chemical (or electrochemical) reaction. 

Model A 

The equations involved in this model can be obtained from the work of 
Bockris et al. (ref. 9). They investigated the electrochemical mechanism for 
deposition of zinc from alkaline solution. Four reactions were found to be 


involved. 

Zn(0H)= I Zn(0H)3 + 0H~ (1) 

RDS 

Zn(0H)“ + e" J Zn(0H)~ + 0H~ (2) 

Zn(0H)2 ^ ZnOH + OH" (3) 

ZnOH + e" ^ Zn° + OH" (4) 


Reaction (2) is the Rate Determining Step (RDS). Because reactions (1) and 
(3) are nonelectrochemical and in equilibrium, the mechanism may be abbre- 
viated for our purposes. Reactions (1) and (2) may be combined into reaction 
(1-2) and reactions (3) and (4) may be combined into reaction (3-4). 

RDS 

Zn(0H)J + e" Zn(OH)" + 2 OH" ' (1-2) 

Zn(OH)" + e" ^ Zn° + 2 OH" 

Looking at the reactions on the basal and non-basal planes individually allows 
one to generate four Butler-Volmer type equations. Such equations now involve 
the area fraction of basal and non-basal planes. To simplify the notation and 
manipulation of these equations, we express them as follows: 
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Ni = N2 (1 - Xn v2)/vNA 
N3 = N4 (Xn - v 2 )/vNB 
Bi = B2 (1 - Xb v2)/vBA 
B3 = B4 (Xb - V2)/VBB 


( 5 ) 

( 6 ) 

( 7 ) 

( 8 ) 


The Ni, N3, Bi, and B3 are currents per unit area; the N's refer to 
the non-basal planes while the B's refer to the basal plane; the subscript 1 
refers to reactions ( 1 - 2 ) and the subscript 3 refers to reactions ( 3 - 4 ). The 
N2, N4, B2, and B4 are the products of the corresponding exchange cur- 
rent and the area fraction of the plane involved. (One should note that the 
areas include only those in the plating region and not regions which do not 
"feel" the plating overvoltage.) Because equation ( 1 - 2 ) is the rate deter- 
mining step, either B2 » B4 or N2 >> N4. The symbol V is equal to 
exp (Fv/ 2 RT) where v is the overvoltage and F is the Faraday. The NA, NB, BA, 

and BB are twice the corresponding symmetry factors. The X^ and Xb are 
the ratios of the concentration of the intermediate species Zn( 0 H )2 at an 
overvoltage v compared with its concentration at v = 0 . The concentration of 
the other species at the surface such as Zn(0H)4 do not appear explicitly in 
equations ( 5 ) to (8) but are implicitly incorporated into the exchange cur- 
rents. The surface concentrations of such species are assumed to be constant, 
independent of voltage. This is an oversimplification because concentration 
polarization does occur. But this simplification does not invalidate the 
results of this paper. _ 

The concentration of the intermediate species Zn(0H)2 on one plane 
is dependent on its value on the other plane because of possible diffusion 
from the plane with the higher concentration to the plane with the lower 
concentration. 


where D2 is the rate constant which is dependent on the areas of each 
individual plane. If the facets of the individual planes are small, D2 will 
be relatively large, allowing rapid diffusion. Again, for simplicity we con- 
sider D2 to be constant for any one particular deposit. Equation ( 10 ) is 
written as if diffusion proceeded from the N planes to the B planes. If 
Xb > Xfj and the opposite occurs, will be negative. D2 must be 
positive. The mathematical details of the model A are given in appendix A. 

The results of our calculations are summarized in figures 4 to 7 . In 
these figures R (N3/N5; plating on non-basal plane divided by plating on 
basal plane) is plotted against the cathodic overvoltages. Each figure (4 to 
7 ) represents a different set of basic constants N?, N4, B2, and B4 
while each subfigure (b) to (j) shows the effect of varying one of the nine 
rate constants. The plots show that the value of R can vary considerably with 


Zn(0H)2 (one plane) » Zn(0H)2 (other plane) 


( 9 ) 


Mathematically the diffusion rate Di can be expressed as 



( 10 ) 
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voltage if an exponential constant of the rate controlling step is varied (See 
figs. 4(d), 4(i), 5(d), and 5(i)). The value of R can even be varied from 
above unity to below unity as a function of voltage if N 2 > B 4 > N 4 > 62 
(See fig. 7). This indicates that at low voltage, plating on the basal 
plane can dominate while at high voltages, plating on the non-basal planes 
can dominate. Moreover, in some cases R can even reach a maximum. Finally 
and most importantly, R cannot have a negative value in model A as implied 
in figures 4 to 7, and as shown algebraicly in appendix A. Model A alone 
would be acceptable if the mechanism did not require dissolution from the N 
planes at low cathodic voltages. 


Other Models 

Other possible models are easily devised by incorporating an auxiliary 
reaction into model A for each. Such reactions may involve disproportiona- 
tion, simple corrosion, electrochemical corrosion, or catalytic redox. Some 
of these reactions are very likely to occur such as those involving electro- 
chemical corrosion (ref. 10); some are not too likely. The 13 we tried are 
listed below. The equation given for a particular model is generally more 
discriptive than the model's name. 

Model B Disproportionation 
2 Zn( 0 H )2 I Zn° + Zn(0H)= 

Model C Nonelectrochemical oxidation 
Zn( 0 H )2 + H^O + OH" Zn( 0 H )4 + H (adsorbed) 

Model D Nonelectrochemical corrosion 
Zn° + H^O + OH" ^ Zn(OH)" + H (adsorbed) 

Model E Coupled chemical and electrochemical oxidation 

Zn(0H)“ + 2 H 2 O + e" Zn(0H)J + 2H (adsorbed) (14) 

Model F Coupled chemical and electrochemical corrosion 
Zn° + 2 H^O + e" ^ Zn( 0 H )2 + 2H (adsorbed) 

Model G Coupled plating and hydrogen formation 
Zn(0H)~ + H^O + 2e" J Zn° + 3 OH" + H (adsorbed) 

Model H Coupled reduction and hydrogen formation 
Zn(0H)J + H 2 O + 2e" ^ Zn( 0 H )2 + 3 OH" + H (adsorbed) 


(15) 

(16) 
(17) 


( 11 ) 

( 12 ) 

(13) 
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Model I Catalytic plating 

2Zn{0H)~ + e~ Zn(OH)" + Zn° + 2 OH" (18) 

Model J Catalytic reduction 

Zn(0H)4 + Zn(0H)“ + e" ^ 2Zn(OH)2 + 2 OH" (19) 

Model K Multiple nonelectrochemical reduction 

3H (adsorbed) + 2Zn(0H)^ ^ Zn° + Zn(OH)" + 3 OH" + OH^O (20) 

Model L Multiple chemical plus electrochemical reduction 
2H (adsorbed) + 2Zn(0H)= + e" Zn° + Zn(0H)~ + 4 OH" + 2H^0 (21) 

Model M Multiple electrochemical plus chemical reduction 
H (adsorbed) + 2Zn(0H)4 + 2e" ^ Zn° + Zn(0H)2 + 5 OH" + H 2 O (22) 

Model N Multiple electrochemical reaction 

2Zn(0H)^ + 3e" ^ Zn° + Zn(0H)~ + 6 OH" (23) 


The mathematical treatments and other details of these models are given in 
appendix B. The mathematical treatments were aimed at finding sets of con- 
stants which would result in a zero plating current on the N planes at some 
cathodic voltage (implying positive values at higher voltages and negative 
values at lower voltages). The Butler-Volmer equations are given in table I 
and the important computer steps are listed in tables II, III, and IV. A wide 
variety of constants were tried (tables V and VI), but for most of the models, 
a zero value for the non-basal plating current required a negative value for 
either a rate constant (exchange current) or a concentration. Thus, models C, 
D, E, F, G, H, I, J, and N were found to be unacceptable. However, for cer- 
tain positive values of constants for models B, K, L, and M, the non-basal 
plating current was zero at cathodic voltages, and the non-basal plating cur- 
rent was negative at lower cathodic voltages. In table VII, the sets and 
subsets for which this is true is indicated by a positive value for Rj. 

Thus, it is possible to have dissolution for certain planes in the cathodic 
region. 

Once one knows that some models will work, it is interesting to examine 
the dependence of the various currents and other variables on voltage. The 
pertinent program steps are indicated in table VIII. We will look at one set 
of conditions for each workable model. A graphic survey of all the conditions 
(sets of constants) for which Np (plating current on the non-basal planes) = 0 
is beyond the scope of this report. 

Model B . - This model involves disproportionation, equation (11), along 
with the fundamental equations (model A). It was found to result in a feasi- 
ble model only when disproportionation is very rapid on the N planes and non- 
existent on the B plane; mathematically, X|\| » 1 (Xj^ = 1) and B 5 = 0. 
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In more practical terms this implies the N planes have a much higher catalytic 
activity than the B plane. 

Figure 8 shows the dependence of various currents on overvoltage for this 
model for one particular set of constants (set 2, subset 0, D 2 = 13.86). 

In this and subsequent figures solid lines indicate currents in the plating 
direction and dashed lines indicate currents in the dissolving direction. 

In figure 8 (a), the plating current Bp on the basal plane is well behaved: 
in the cathodic region plating occurs, and in the anodic region dissolution 
occurs. However, the plating current on the non-basal planes Np exhibits 
the type of anomolous behavior that we have been seeking: at low cathodic 

voltages, dissolution occurs. At higher cathodic voltages and at anodic volt- 
ages, plating occurs on the N planes. Plating on the N planes at anodic volt- 
ages would be very difficult to observe. 

It is easier to examine higher current - voltage behavior via the Tafel 
type plot in figure 8 (b). In this figure, one sees that at sufficiently high 
cathodic voltages, the plating on the N planes exceeds that on the B plane due 
to an inflection of the slope of Bp in the 100 to 200 mV region. This sat- 
isfies the condition needed for dendrite formation. From figure 8 (b), one 
also sees that at high anodic voltages the anomolous plating on the N planes 
in the anodic region finally changes to dissolution. 

Figure 8 (c) shows the overall Tafel plot that would be observed for the 
system; the log of the total current is plotted against the voltage. This 
illustrates that the anomolous behavior of individual planes would not show up 
in a Tafel experiment. Moreover, the Tafel slopes (fig. 8 (c)) are in agree- 
ment with experimental values (ref. 8 ). At this time, we point out that the 
shape of the curves depend only on the ratio of rate constants and not on 
their absolute values. The absolute values affect the position of the curves 
on the current scale; that is they affect the value of the exchange current, 
ip. Thus, any model can be adjusted so as to give an iQ in agreement with 
the experimental data. 

Finally, figure 8 (d) shows the dependence on voltage of the individual 
non-basal plane current components: N 3 , and N 5 . The currents 

and N 3 exhibit conventional Tafel behavior at all voltages. It is N 5 that 
is responsible for all the anomolous behavior; N 5 is the cause of dis- 
solution in the cathodic region and plating in the anodic region. 

Models K, L, and M. - The chemistries of these three models are quite 
similar as seen from a comparison of equations (21) to (23). They differ only 
in the relative amounts of hydrogen and electron involvement in their auxil- 
iary reaction, having in common all the fundamental plating reactions 
(model A). Model K is unique with respect to the other two in that its 
auxiliary reaction (eq. ( 21 )) is not electrochemical. 

Plating currents versus voltage curves are presented in figures 9, 10, 
and 11, respectively, for models K, L, and M. The sets of constants used are 
indicated in their respective table legends. The behavior shown for models K, 
L, and M is very similar to the type exhibited in figures 8 (a) and (b) for 
model B. Current for the basal plane is again well behaved: cathodic in the 

cathodic region and anodic in the anodic region. The current for the non- 
basal plane is anodic in the low cathodic voltage region and cathodic in the 
low anodic voltage region. In higher cathodic voltage regions, the N current 
becomes cathodic. Model K at high anodic regions shows anodic behavior for 
the N planes which makes model K in complete qualitative agreement with model 
B as far as current voltage dependence is concerned. Perhaps models L and M 
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were not taken to a high enough anodic voltage for the N planes to exhibit 
anodic behavior in the anodic region. The very close similarity between 
models B and K could be attributable to the fact that both these models have 
nonelectrochemical auxiliary reactions, equations (11) and (20). 

Although not shown, models K, L, and M exhibit conventional Tafel curves 
as does model B. And, for each of these models, the anomolous non-basal plane 
behavior is due to the N5 current being in the opposite direction from 
N3. Thus, at this point there is no grounds for a selection of a "best" 
model from among B, K, L, and M. 

Effect of hydrogen . - There is one major difference between model B and 
models K, L, and M. Model B does not involve the presence of hydrogen 
(eq. (11)); the other three do (see eqs. (21) to (23)). Therefore, we 
examined the hydrogen activity for models K, L, and M as a function of volt- 
age; table VIII outlines the procedure. Figure 12 shows the variety of 
dependencies that can be obtained. Hydrogen activity can increase with 
cathodic voltage and decrease with anodic voltage as in the case of model K. 
The activity can increase with both cathodic and anodic voltage; see condi- 
tions K'. Finally for models L and M, the hydrogen activity can increase with 
anodic voltage and decrease with cathodic voltage. If this is representative 
behavior of models L and M, for all combinations of constants, these models 
are not acceptable because hydrogen pressure (activity) is known to increase 
in the cathodic region of the alkaline zinc electrode. In contrast, model K 
is acceptable in that it predicts hydrogen formation at high cathodic voltages 
(dendritic region) as is generally observed. It has been reported (ref. 11) 
that the alkaline zinc electrode also evolves hydrogen at high anodic volt- 
ages. If this be the case model K is still acceptable because conditions K' 
predict the increased hydrogen pressure in both the anodic and cathodic 
regions. 

Any of the models which involve hydrogen activity in their reactions will 
evolve hydrogen gas when the value of H corresponds to a pressure above atmos- 
pheric. This condition sets an upper limit to the value of H. For figure 13, 
the assumption was made of a fixed maximum for H for model K with conditions 
K' (for this condition hydrogen is produced in both the anodic and cathodic 
directions). The effect is that in the anodic region, Np switches from 
plating to dissolving at a lower voltage. Of more importance is the fact that 
a negligible effect occurs in the cathodic region; thus, ignoring the fact H 
may reach a maximum value does not affect any of our previous conclusions. 

CONCLUDING REMARKS 

From an examination of the appearance and crystal orientation of the vari- 
ous zinc deposits, we have postulated the following mechanism as an explana- 
tion for the variety of morphologies. The dendrite morphology produced at 
high current densities is a result of more rapid deposition on the non-basal 
crystallographic planes (ITOO, lOTO, and OlTO) than on the basal plane. 

The mossy morphology produced at low current densities is due to actual dis- 
solution of zinc from the non-basal planes at low cathodic voltages. The 
boulder modification which forms at intermediate current densities results 
from comparable plating on both the basal and non-basal planes. 

Fourteen electrochemical models were explored to account for such behav- 
ior. In all of these, the exchange currents, rate constants, etc., were 
allowed to have a variety of values and they were assumed to be different for 
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the basal and non-basal planes. The currents versus voltage were calculated 
via computer programs. A model consisting only of fundamental plating reac- 
tions can explain different rates on different planes but not dissolution from 
a plane in the cathodic voltage region. Each of the other models consists of 
the fundamental plating reactions plus an additional reaction. Only two of 
the models explored can account for the proposed mechanism. The one model 
involves very rapid disproportion of the Zn(0H)2 intermediate to form 
zincate and zinc on the non-basal planes and no disproportionation on the 
basal plane. The second workable model involves a redox reaction between the 
zinc-zincate and hydrogen-water species. The second model can thus be used to 
explain the evolution of hydrogen at high cathodic voltages. It can also 
explain the hydrogen evolution at anodic voltages reported (ref. 11) for the 
alkaline zinc system. The two workable models give overall Tafel data in 
agreement with experimental work (ref. 9). 

The general mechanism and models may be applicable to other metal systems 
with multiple morphologies. The metal involved must be at least divalent. 
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APPENDIX A - MODEL A 


We again re-emphasize all the models assume steady state conditions. 

Model A involves only diffusion of the intermediate species, reaction ( 9 ), and 
the electrochemical zinc plating, reactions ( 1 - 2 ) and ( 3 - 4 ). Thus, equations 
( 5 ) to (8) and equation ( 10 ) are involved. Equation ( 10 ) requires a little 
more discussion than given in the text proper. 


Di = D2 m - Xb) (10) 

For ease of manipulation, the D^ is expressed in units of current per unit 
area to agree with the units in equations ( 5 ) to (8). D£ has to be ex- 
pressed in the same units. The D2 comprises the diffusion coefficient, the 
length of border between basal and non-basal surface planes, and the surface 
areas of both sets of surface planes. In this paper we will not express D2 
in more detailed terms because of the complexities that would arise. 

The important steps in the computer program are as follows: 


Assign values to the constants (Al) 

For V = -10 to -390 Step- 20 : Rem V = millivolts; 

negative voltages are cathodic (A 2 ) 

V = exp (+0.01919V) (A 3 ) 

87 = B4/VBB + 02 v(2-BA) (A4) 

Bs = 82 /vBA + B4 v( 2 -BB) (A 5 ) 

N7 = N4/V'^B + 1^2 v( 2 -NA) (A6) 

Ns = N2/V^^ + N4 (A 7 ) 

Di = D2 (N8/N7 - Bs/B7)/(1 + D2/N7 + D2/B7) (A8) 

Xn = (N 0 - Di)/N7 (A 9 ) 

Xb = (Bs + Di)/B 7 (AlO) 

Equations ( 5 ) to (8) (All) 

R = N3/B3 (A 12 ) 

Print values (A 13 ) 


To study model A, we have used four basic sets of constants that were in 
keeping with the Bockris electrochemical mechanism. Values of R were cal- 
culated as a function of overvoltage. Then, the effect of changing the value 
of one constant at a time was recorded. The plots of R versus v thus obtained 
are given in figures 4 to 7 . 

Some of the conclusions are given in the text proper. More details are 
given here. Some of the conclusions are obvious without even reference to 
these figures; some are not. When the N parameters (N2, N4, NA, NB) are 
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equal to the respective B parameters, R is unity independent of voltage 
(fig. 4(a)); and there is no dependence on D 2 (fig. 4(f)). If changing an 
N parameter has an effect, then changing the respective B parameter the same 
amount has about the same magnitude of effect but in the opposite direction. 
Compare figure 4(b) with 4(g); 4(d) with 4(i); 5(d) with 5(i); etc. Changing 
a constant for the fast step (one with a relatively high constant) has little 
or no effect on R, or on its dependence on v (see figs. 4(c), 4(e), 4(h), 

4(j), 5(c), 5(e), 5(h), 5(j), etc.). Changing the ratio between the constant 
for the slow reaction on one plane with respect to that on the other plane 
changes the value of R; compare figure 4(b) with 4(a), 4(g) with 4(a), 5(b) 
with 5(a), 5(g) with 5(a), 6 (b) with 6 (a), etc. For these cases in figures 4 
and 5, R becomes equal to Np/B?; for these cases in figures 6 and 7, R 
becomes equal to N 2 /B 2 at the higher voltages. Appreciable changes in R 
with voltage arise from two conditions: ( 1 ) variations of an exponential con- 

stant of the rate determining step (NA and BA (figs. 4(d), 4(i), 5(d), and 
5(i)); and (2) the condition of having N 2 > N 4 while B 4 > B 2 (figs. 6 
and 7). Note that in figure 7 when > N 2 > Bo > N^, the values of R 
do cross the value of unity at some cathodic voltage. The curves in figure 7 
indicate that at low voltages, plating on the basal planes can dominate and at 
high voltages, plating on the non-basal planes can dominate. If both an 
exponential constant varies from unity and N 2 > N 4 while B 4 > Bp, then 
the value of can reach a maximum at some voltage. (See figs. 6 (i) and 7(i).) 

In order for R to have a negative value, it must equal zero at some 
cathodic voltage. Thus, 


N 3 = 0 

Substituting this into equation ( 6 ) 


(A14) 


Xn . V2 

Due to conservation of Zn( 0 H )2 on each plane 

Di = Ni 


(A16) 


(A15) 


and 

Bi + Di = B 3 

Substituting equations (5), (7), ( 8 ), and (10) into these equations 
Xb = V2 + (N 2 (vA - 1))/(D2 VNA) 


(A18) 


(A17) 


and 


Xb = (Bb + D2 v2)/(B7 + D 2 ) 

The symbols B 7 and Bg are defined in equations (A4) and (A5). 
Combining these: 


(A19) 


N 2 /D 2 = ((Bb + D 2 v2)/(B7 + D 2 ) - V2) vNA/(v4 _ 1) 


(A20) 
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(A21) 


Because N 2 /D 2 > 0 and _ 1 < 0 , then 
(Bs + D2 v2)/(B7 + D 2 ) < V2 

Simplifying and substituting the values of B 7 and Bs from equations (A4) 
and (A5) 


V2 > (B 4 /VBB + B 2 v( 2 -BA))/(B 2 /vBA + B 4 v(2-BB)) (a22) 

Additional simplification gives 

> 1 (A23) 

But V < 1 in the cathodic region, so that No ^ 0 in the cathodic region 
and thus for model A, there cannot be dissolution of metal on the non-basal 
plane in the cathodic region. 
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APPENDIX B - OTHER MODELS 


Each model involves the reactions for model A; thus, for each model equa- 
tions ( 1 ) to ( 10 ) are applicable. In addition, each model includes one of the 
reactions, equations ( 11 ) to ( 23 ). As a result, each model also has two addi- 
tional mathematical (Butler-Volmer) relationships as given in table I: one 

for the basal planes and one for the non-basal planes. The N5 and B5 are 
the respective rates on the non-basal and basal planes. These are given in 
current per unit area for the sake of agreement with other equations used. 

The Np is the total plating current on the N plane. The N5 and B5 are 
the corresponding rate constants. The NC and BC are equal twice the corre- 
sponding symmetry factors. The same set of symbols is used in each model to 
avoid complexities in symbol notation. The symbol H refers to the ratio of 
hydrogen activity at voltage v to its value at v = 0 . The diffusion of hydro- 
gen absorbed on the electrode surfaces was assumed to be much more rapid than 
that for the Zn(0H)2 species so that the value of H was taken to be the 

same on both sets of planes. The expressions for N^, 85, and Np were 

obtained by examining the respective chemical reaction, equations (11) to 
( 23 ), as well as equations ( 1 - 2 ) and ( 3 - 4 ). 

From the expressions in table I, equations ( 5 ) to (8) and ( 10 ), and the 

condition of steady state; computer programs were set up to calculate a ratio 

Rj between two selected rate constants under the condition Np = 0 . The 
important steps in these programs are listed in table II for D2 = 0 . In 
general permitting D2 > 0 would allow Xj^j to become closer to Xg in value 
and thus would nullify in part the effect of the auxiliary equations intro- 
duced into these models. 

Only for models B, I, J, and N did the introduction of diffusion appear to 
be a possible way of improving the model, that is allowing it to exhibit 
Np = 0 for some set of rate constants. We thus set up for models B, I, J, 
and N, computer programs for which D2 was very large so that X|\j = Xg: 
the essence of these programs are listed in table III. 

For model B still another possibility exists because the auxiliary equa- 
tion neither involves an electrode reaction nor the presence of the "foreign" 
species, hydrogen. Therefore, in model B, B5 does not have to equal -N5 
(a necessity when H is produced at one surface used up at the other). We let 
65 = 0 and let Ns have a very large value by setting X^ = 1, meaning 
Zn(0H)2 was in equilibrium with Zn° and Zn(0H)4. The important computer 
steps are in table IV. 

One must remember that Np = 0 is not the only requirement for a feasible 
model. It is also necessary for (aNp/av)^]^ to be positive where Np = 0 . 

To determine this we made use of the following equation, all derivatives being 
evaluated at Np = 0 . 

(aNp/av)pi = - ( aRj/ 3v)Np/ ( 3 Ri/ aNp)v 

Into each computer program, we inserted a subroutine involving the presence of 
a slight amount of plating on the non-basal plane, Np, and thus determined 
whether (aRj^/aNp)^ was positive or negative. The dependence of R^ on 
V from our original program gives us the sign of (aR]^/av)Np. In addition, 
we added another subroutine to ascertain if Xj^, Xg, and H (hydrogen con- 
centration in appropriate models) had positive values. 
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It is readily seen from the final equation in some of the programs (tables 
II, III, and IV) that the ratio Ri between the rate constants determined is 
always negative regardless of the values of other constants. Thus, for these 
models, Np can not be zero in the cathodic region. These models have been 
indicated by asterisks in tables II, and III; the models are models I and J 
with = 0 and model I and J with = Xg. The programs of all the 
models were run using nine different sets of constants as shown in table V. 

By close examination of the equations one can determine that the actual value 
of the constants is not important in the calculation, only their relative 
values. For each set, each program was rerun with the change of one constant 
with respect to the original set. The changes are indicated in table VI. 

By this method we found that for models B, C, D, E, F, G, H, I, 0, and N with 
D 2 = 0 and for models B, I, and N with Xf^ = Xg, it was not possible to 
have Np = 0 in the cathodic region. 

However, for some models and certain values of the constants, it was 
indeed found possible to have Np = 0 in the cathodic region: these included 

models K, L, and M with Do = 0; ahd model B with X|^ = 1 and B 5 = 0. The 
sets and subsets for whicn this is true (Np = 0) are indicated in table VII 
by giving the value of Rj for the condition v = -10 mV (10 mV cathodic). 

The programs showed that Np is negative between 0 and -10 mV for the value 
of R^ given in table VII. Some of the sets and subsets are redundant for a 
particular model because (as mentioned earlier) only the ratio of constants 
and not their actual values are important in finding the condition of Np = 0. 
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TABLE I. - EQUATIONS FOR MODELS B TO N 


Model 

^5 

Expressions for 

Np 

B 

Ng (X2 - 1) 

Bg (X^ - 1) 

Ns^Ng 

C 

"6 (Xn - 

Be (Xb - «) 

^3 

D 

Ne (1 - H X^) 

Bg (1 - H Xg) 

^ 3-^5 

E 

Ng (X^ - V^)/V'^^ 

Bg (Xg - V^)/V®^ 

N3 

F 

Ng (1 - 

Bg (1 - Xg V^)/V^^ 

'3 -'5 

G 

Ng (X^ - H 

Bg (Xg - H 

N3- Ng 

H 

Ng (1 - H X|^ 

Bg (1 - HXg 

^3 

I 

Vn 

Vb (’^b ■ 


0 

¥n -^N 

BgXg (1 -Xg V^)/V^^ 

'3 

K 

h - ^n) 

Bg (hB - Xg) 

'3 "'5 

L 

Ng (H^ - V^)/V'^^ 

Bg (H^ - Xg V^)/V®^ 

N3 ^ Ng 

M 

Ng (H - X^ V^)/V^'^‘^ 

Bg (H - Xg 

N3 ^ Ng 

N 

Ng (1 - X|^ 

Bg (1 - Xg 

N3 - Ng 
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TABLE II. - PROGRAMS WHEN 02=0 


[Those models marked with * can be seen to be invalid because 
number. ] 

Model B 

Equations (Al) to (A7) 

^N " - N^/V^^) 

Xg = (-B7 + sqr (B^ + 8 Bg (Bg + 2 Bg))/4 Bg 
R^ = (Ng/N^) = (X^ - v2)/(v'^^(l - X^ 2)) 

Model C 

Equations (Al) to (A7) 



= Ne " ^6 " 

B. NgX^* (Bg Bg/(B^ * Bg)) 

H = B/A 

R^ = (Ng/Ng) = (1 - X^)/(V'^'^ (X^ - H)) 

Model D 

Equations (Al) to (A7) 



A = Ng Bg X, 

B . Ng B 7 * Bg(Bg - Ng) 

C - - B^(Ng * Bg) 

H . (-B + sqr (B^ - 4 AC))/2A 

«i = (VV ' (’'n - ■ '‘n "> 

, 

I 

; 

I 


= neg. 


I 
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TABLE II. - CONTINUED. 


Model E 


Equations (Al) to (A7) 



A = (Ng/v”^^ + Bg/V^^ - B^)/(B^ + Bg) 

B = Ng X^/V'*'^ - Bj; * Bj) 

H = sqr (B/A) 

Rl = (Ng/N^) = (1 - V^) V^) 


Model F 


Equations (Al) to (A7) 

X - V~^ 

N ■ 

^ = Nfi B. X. 

6 6 N 

B = Ng B 7 X^ + Bg (V^ " (Bg - Ng/V"^^) 

C = - B 7 (Bg/V®^ + Ng/V*^^) 

H = sqr ((-B + sqr (B^ - 4 A C)) / 2 A) 

Rj = (Ng/N^) = (X^ - V^) - X^hV) 


Model G 


Equations (Al) to (A7) 



W = 

A = W^(Ng/w'^^ + Bg/W®^ - B) 

B = NgX^/W*^^ + BgBg/(B7wBC + Bg) 

H = B/A 

Rl = (N5/N4) = (V^ - X^) W'^^/(V'^^(X^ - HW^)) 
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Table II. - CONTINUED. 


Model H 

Equations (Al) to (A7) 

W = 

A = Xn Ng Bg - BC - NC) 

B = Ng B? + Bg (Bg + Ng/w'^^) 

C = - By (Ng/w”^^ + Bg/W^^) 

H = (-B + sqr (B^ - 4AC))/2A 

Rj = (Ng/N^) = (X^ - 1) (1 - X^ W^) 


Model I 


Equations (Al) to (A7) 

X - V“^ 

^N ■ 

*R^ = (Ng/N^) = (V^ - X|^) v'^^-^®/(X^ - V^) = - 


Model J 


Equations (Al) to (A7) 

2 

*Rj = (Ng/N^) = (X^ - 1) - X^ V^) = - 
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TABLE II. - CONCLUDED. 


Models K, L, M, and N 

Equations (Al) to (A7) 

If K, then Q = 3; If L, then Q = 2: If M, then Q = 1: If N, then Q 

W = Q 

Xfj = (2N4 + N2 /v'^^)/{2N4/V^® + ~ 

B = Ng + Bg Bg “ ^^/(By + Bg 

A = Ng/W^^ + B/W^^ - Bg ■ ^®^/(B7 + Bg ~ 

If Q > 0, then H = (B/A)^^^ 

If 0 = 0, then H = 1 

Rj = (Ng/N^) = (X^ + V^)/(V'^^ (H^/W^^ - X^ W'^^)) 


TABLE III. - PROGRAMS WHEN = Xg 
Model B 

Equations (Al) to (A7) 

A = 2 Bg + Ng 
B = B 7 + 

C = - (Bg + + A) 

X^ = (-B + sqr (B^ - 4AC))/2A 

^B = 

Rl = (Ng/N^) = (X^ - V^)/(V'^^ (1 - X^ 2)) 


Model I 

Equations (Al) to (A7) 

A = Bg/V®^ 

B = B^/V^® - Bg + N 2 + B^ 

C = B^ - N2/v'^^ - B2/V^^ 

X^ = (-B + sqr (B^ - 4AC))/2A 

'b = 'n 

*Rj = (Ng/N4)=(v2 - X^) - V^)) = - 


Model J 


Equations (Al) to (A7) 

X - 

^B " 

*«!= (N6/N2)=-(v'^^/N2)(Bg/V®‘^+B2/(X^V®^)+N2/(X^v'^^)) 
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TABLE III. - CONCLUDED. 


Model N 

Equations (Al) to (A7) 

W = 

A = + BgW^ - BC - BA + - NC + 

B = B^/V^ - BB + Bg/W®*^ + B2/V^^ + 2Ng/v’^‘^ + N2/v'^^ 
= B/A 

^B " 

*R^ = (Ng/N^) = (V^ - X^) (1 - W^)) 


- NA 
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TABLE IV. - PROGRAM WHEN Bg = 0 AND 
Model B 

Equations (Al) to (A7) 

Xn = 1 

N^ = N 2 (1 - V^)/V^^ 

N3 = N^ (1 - V^)/V'^^ 

= Ni - N3 

Xg = (Bg - D^)/B7 

= Dj/(1 - Xg) 
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TABLE V. - SETS OF CONSTANTS USED 


Set 

B2 

B4 

N2 

N4 

Value of all other constants 

1 

1 

100 

1 

100 

1 

2 

1 

100 

.1 

10 

1 

3 

1 

100 

10 

.1 

1 

4 

1 

100 

100 

1 

1 

5 

100 

1 

1 

100 

1 

6 

100 

1 

.1 

10 

1 

7 

1 

100 

1000 

10 

1 

8 

100 

1 

100 

1 

1 

9 

100 

1 

10 

1000 

1 
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TABLE VI. - VARIATIONS OF CONSTANTS 


Subset 

Value 

0 

As indicated in table I 

1 

N2 = 10 times value of N^ in set 

2 

N4 = 10 times value of N^ in set 

3 

NA = 0.8 

4 

NB = 0.8 

5 

^2 = 10 times value of B^ in set 

6 

8^ = 10 times value of in set 

7 

BA = 0.8 

8 

BB = 0.8 

9 

NC = 0.8 

10 

Bg . 10 

11 

BC = 0.8 
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TABLE VII. - VALUES OF Ri FOR WHICH Np = 0 AT v = 10 mV 


[- indicates negative values] 


Model B (B5 = 0 and = 1) 

^Subset 

Set^ 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

1 

7 . 

3 

4 

5 to 9 

13.86 

13.86 

15.25 

15.28 

13.85 

6719 

13.29 

6719 

13.29 

13.85 

16.36 

16.36 

136.4 

12.38 

136.3 

12.38 

13.85 

13.85 

13.93 

13.93 

^13.86 

13.86 

13.86 

13.86 

13.86 

13.86 

Model K 

1,2 

3,4,7 

5,6,9 

8 

?.450 

3.062 

2.450 

210.6 

2.455 

2.450 

2.358 

2.603 

3.032 

2.426 

23.79 

2.449 

2.451 

2.451 

1.352 

2.450 

Model L 

1,2 
3,4,7 
5, 6, 8, 9 

2.966 

3.707 

2.966 

2.972 

2.966 

2.855 

3.152 

3.673 

2.939 

2.966 

2.967 

.07653 

2.328 

1.636 

3.024 

Model M 

1,2 

3,4,7 

5,6,9 

8 

3.592 

4.489 

3.591 

308.6 

3.598 

3.592 

3.456 

3.817 

4.450 

3.560 

34.86 

3.591 

3.593 

.04574 

2.319 

6.268 

1.980 

3.735 


TABLE VIII. - STEPS FOR CALCULATING CURRENTS 
VALUES AS A FUNCTION OF VOLTAGE 
Model B (Xj^i = 1; B5 = 0) 

Assign values to constants 

Assign value to D2 which makes Np = 0 at v = -10 mv 
For V = -30 to +30 step 2 [or for -300 to +300 step 20] 
Equations (A3) to (All) 

Appropriate expressions for 85 and Nfi from table I 
Print Results 


Model K, L, and M (D2 = 0) 

Assign values to constants 

Assign value to Rj which makes Np = 0 at v = -10 mv 
Let N4 = R^ N5 

Let N2 = N4 times original value of N2/N4 

For V = -30 to +30 step 2 [or for -300 to +300 step 20] 

Equations (A3) to (All) 

Appropriate expressions for B5 and Ns from table I 
Print Results 
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Figure 1. - SEM view of dendritic zinc deposit. 




(a) SEM VIEW OF MOSSY ZINC DEPOSIT. 



0.2 pm 


(b) TEM VIEW OF MOSSY ZINC SINGLE CRYSTAL 
Figure 2. 





(C) ELECTRON DIFFRACTION PATTERN OF ZINC MOSS CRYSTAL OF FIGURE 2(b). 
Figure Z - Concluded. 



Figure 3. - 


SEM view of zinc boulder deposit. 
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Figure 4. - Ratio (R) of plating on nonbasal plane to basal plane as a function of voltage: Model A O locates R - 1; A locates -200 mV; unless 
otherwise specified; 63 ■ 1; B 4 - 100; BA • 1; BB - 1; N 2 ■ 1; N 4 - 100; NA • 1; NB - 1; and D 2 • 0). 


> 

> 

> 

> C 

^ ( 

) c 

) ( 




(a) AS SPECIFIED. 

(b) B 2 - 0 . 1 . 

(c) BA - la 

(d) BA - a 9. 

(e) BB ■ 0.9. 

> 

> 

> 

> 

> 

i { 

■ <1 

» ( 

) ( 

) ^ 


(f) D 2 ■ 100. (g) N 2 - 0. 01. (h) N 4 • 1. (I) NA - a 9. (j) NB • a 9. 

Figure 5. - Ratio (R) of plating on nonbasal plane to basal plane as a function of voltage: Model A Olocates R * 1; A locates -200 mV; unless 

otherwise specified; B 2 -1; 84 ■ 100; BA - 1; BB • 1; N 2 - a 1; N 4 ■ 10; NA ■ 1; NB ■ 1; and D 2 ■ 0). 
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(i) NA - a 9. 


(e) 88 - 0 . 9. 
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Figure 6 . - Ratio (R) of plating on nonbasal plane to basal plane as a function of voltage-, Model A Olocates R • 1; A locates -200 mV; unless 
otherwise specified; 83 • 1; 84 ■ 100; 8 A - 1; 88 - 1; N 2 - 100; N 4 - 1; NA - 1; N 8 - 1; and D 2 • 0). 
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Figure 7. - Ratio (R) of plating on nonbasal plane to basal plane as function of voltage-. Model A Olocates R - 1; A locates -200 mV; unless 
otherwise specified: 82 -1; 84 ■ 100; 8 A • 1; 88-1; N 2 - 3. 0*. N 4 • a 03; NA ■ 1; N 8 ■ 1; D 2 ■ 0). 


PLATING CURRENT TOTAL CURRENT PLATING CURRENT 




(a) PLATING CURRENT AT LOW (b) PLATING CURRENT AT HIGH 

VOLTAGE. VOLTAGE. 



(c) ALGEBRAIC TOTAL OF ALL (d) COMPONENTS OF THE NON- 
CURRENTS. BASAL CURRENTS. 

Figure 8. - Current vs. voltage curves for Model B. (Solid line - plating direc- 
tion; dashed - dissolving direction; Bo -1. B^ • 100, N 2 • 0. 1, N 4 • 10, B^ - 0, 
N^ - BA - BB - NA • NB • 1, X|^ - 1. D 2 “ 13. 86; rate constant units are 
arbitrary, same as currents.) 




V, mV 

(a) FOR LOW VOLTAGE. (b) FOR HIGH VOLTAGE. 

Figure 9. - Plating current vs. voltage curves for Model K. (Solid line Indicates 
current in the plating direction; dashed line Indicates current in the dissolving 
direction; Bo "1, B 4 - 100, B^ - N^ - 1, BA • BB ■ NA ■ NB - 1, D 2 " 0, 

Rj - 2.450, N 4 - N^/Rj, and N 2 ■ 100 N 4 ; rate constant units are arbitrary, 
but the same as for currents. ) 


PLATING CURRENT. 




(a) FOR LOW VOLTAGE. (b) FOR HIGH VOLTAGE. 

Figure la - Plating current vs. curves for Model L (Solid line indicates current 
in the plating direction; dashed line indicates current in the dissolving direction,- 
B 2 -I. B 4 -IOO, B^- Nj- 1. BA- BB - NA- NB- 1. D^-0, Ri-i% 6 , 

N 4 - N^/Rj, and N 2 • 100 N 4 ; rate constant units are arbitrary, but the same as 
for the currents. ) 



V, mV 



(a) FOR LOW VOLTAGE. (b) FOR HIGH VOLTAGE. 

Figure 11. - Plating current vs. voltage curves for Model M. (Solid line 
indicates current plating in the plating direction; dashed line indicates 
current in the dissolving direction; B 2 ■ 1 . B 4 - 100 , B 5 - N^ • 1 , 

BA - BB • NA - NB • 1. Do • 0. Ri - 3. 592. N 4 • N^/Ri, and N 2 • 100 N 4 ; 
rate constant units are arbitrary, but the same as for the currents.) 



Figure 12, - Hydrogen activity rela- 
tive to its equilibrium value. The 
constants for curves K. L, and M 
are given respectively in figures 9. 
10, and 11. constants are those 
of Model K, set 1, subset 6. 
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Figure 13. - Effect of a maximum hydrogen activity on No for Model K 
(set 1, subset 6). 
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